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Abstract 

The unitary operators U{t), describing the quantum time evolution of systems with 
a time-dependent Hamiltonian, can be constructed in an explicit manner using the 
method of time-dependent invariants. We clarify the role of Lie-algebraic techniques 
in this context and elaborate the theory for SU(2) and SU(1,1). 

We show that the constructions known as Magnus expansion and Wei-Norman 
expansion correspond with different representations of the rotation group. A simpler 
construction is obtained when representing rotations in terms of Euler angles. 

The many applications are postponed to Part II of the paper. 

1 Introduction 

Assume that a time-dependent Hamiltonian H{t) is given. The axioms of quantum mechanics 
require that it is a selfadjoint operator in Hilbert space. The solution of the time dependent 
Schrodinger equation 

= H{t)m (1) 

is formally given by iplt) = U{t)ilj{0) with unitary operators U{t) satisfying U{0) = I and 

t(^^U{t)^U*{t) = H{t). (2) 

The time evolution can be extended to density operators (positive trace-class operators with 
trace 1) by the relation 

pit) = uit)pm*{t). (3) 
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They satisfy von Neumann's equations of motion 

^p(t)=.[p(t),/J(t)]. (4) 

Explicit expressions for U{t) in function of the given H{t) are usuaUy difficult to obtain. 
The mathematical origin of the difficulties is that Hamiltonians H{t) and H{t') at unequal 
times t ^ t' do not necessarily commute with each other. Much progress has been made 
in the case that the time-dependent Hamiltonian can be written as a linear combination of 
constant operators Si, S2, ■ ■ ■ Sn, generating a finite dimensional Lie algebra 

n 

H{t) = Y,hAt)S,. (5) 

i=i 

The generators satisfy the commutation relations 

[Sj, Sm] = i''^^UJjmlSl, (6) 

with structure constants Ujmi- Then the time evolution of the generators can be written as 

U{t)S,U*{t) = J2^mj{t)S^, (7) 



with coefficients Umj{t) solution of 

dt 



^U^j (t) = ^ ^^Ismhl it)Usj (t) . (8) 



Is 



Note that ([7]) is not the time evolution in the Heisenberg picture. The latter is given by 
U*{t)SjU{t) instead of U{t)SjU*{t). 

Using ([7]) one can obtain the general solution of the von Neumann equation (jl]), with 
arbitrary initial conditions of the form 

p{t = 0) = C + J2<^A0)Sj. (9) 
j 

The operator C is an arbitrary operator commuting with all generators Sj. The solution is 

p{t) = c + J2^At)Sj, (10) 

j 

with 

= '^Ujm(t)am{0). (11) 

m 

The problem remains to obtain an explicit expression for the unitary operators U (t). See 
for instance [H [2] and the references quoted there. These authors show that one can write 
U (t) as an ordered product 

U{t) = e^iW^ie^^^W^^ • • • e^'"(*)^", (12) 
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where the complex functions gi{t), ■ ■ ■ gn{t) are solutions of some set of non-linear differential 
equations. In this way the problem is reduced to a well-known but rather difficult classical 
(this is, non-quantum) problem of solving sets of equations. This method has been elaborated 
by many authors, including [3llll[5l[6l[71IHl[9]. A related topic is that of the superposition 
principle for nonlinear equations, see for instance [121 [E] • 

Alternatively, one assumes that a particular solution Ps{t) of von Neumann's equation 
(jl]) is known. This special solution was called a time-dependent invariant in [3]. We use 
the existence of such special solutions as a starting point for constructing unitary operators 
U (t) , based on an idea of \l3l ■ 

The next Section describes in general terms the method of constructing U{t). Section 
3 explains how the assumption, that the Hamiltonian H{t) can be expanded in terms of 
generators of a Lie algebra, helps to apply the method. Subsequent Sections specialise the 
method for SU(2) and SU(1,1). The discussion of the obtained results follows in Section O 
In part II of the paper some applications are considered. 



2 Method 

Assume that there is given a time-dependent Hamiltonian H{t) and a special solution Ps{t) 
of the von Neumann equation (jl]). Can one construct unitary operators U{t) solving ([2])? 
The method described below requires two steps to do so. 

First step Determine unitary operators V{t) such that 

Ps{t) = V{t)ps{0)V*{t). (13) 

If ps(t) has a discrete spectrum, as is the case for a density operator, then this can be 
obtained in principle by diagonalising Ps{t). The unitary V{t) transforms the basis in which 
Ps{0) is diagonal to the basis where ps(t) is diagonal. 

With this V{t) one can construct a Hamiltonian K{t) by 

m = z^v*{t). (14) 

In general, this Hamiltonian differs form the given H{t), in which case V{t) is not a solution 
of ([2]), which is the defining equation of U{t). 

Second step Assume now that ps{0) is self-adjoint, as is the case for a density operator. 
Then generically operators commuting with Ps(0) are functions of Ps(0) in the sense of 
spectral theory. Comparing (fT3|) with ([3]) shows that U(t)V*(t) commutes with ps(t) and 
V*{t)U{t) commutes with Ps(0). Hence, one can expect that there exists a real function ft{x) 
such that 

U{t) = e'^*(^^(*)V(t) = \/(t)e'^*(^^(°)\ (15) 
It remains to determine this function ft{x). Combine (fT5l) with ([2]) to obtain 

^/*(P.(0)) = U*{t) {K{t) - H{t)} U{t) = V*{t) {K{t) - H{t)} V{t). (16) 
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The r.h.s. of this equation is known from the first step of the method. Hence the function 
dft{x)/dt can be determined. To this purpose, the r.h.s. of (ITB]) has to be written as a 
function of Ps{0). Next, /t(x) is obtained by integration. 



3 Lie algebras 

Assume now that H{t) is a hnear combination of generators Si, - ■ ■ , Sn of a Lie algebra, 
with commutation relations as given by (Q, and that the special solution Ps(t) is of the form 
(fTOj) . We first study the time evolution of the coefficients ajit) as elements of a Lie algebra, 
represented in M". 

Introduce the Killing form 



{a\b) = {b\a) = - I ^C^imnOj | I '^^^Inmk 

mn \ j / \ I 



(17) 



The minus sign is needed because of the imaginary factor i in the definition of the 
structure constants. Introduce the Lie bracket a x b = —b x a hj 

(g X b)i = y^^ajbmUjmi- (18) 

jm 

A well-known property of the Killing form, consequence of the Jacobi identity, is that 

(a X b\c) = {a\b x c). (19) 

This will be used often. 

The Hamiltonian K{t) can be written as 

m = Y,msr (20) 



One has 



i^kj{t)am{t)ujjmiSi. (21) 



jm 



This implies 



^ = ^k^it)aiit)ujmij- (22) 



ml 



Introduce the notation a for the vector with components daj/dt. From (122|) then follows 
that a = k X a. 

Let us assume that the Killing form is non-degenerate. This is known to be the case if 
and only if the Lie algebra is semi-simple. 
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Lemma 1 Assume that the Killing form is non-degenerate. Fix a G such that {a\a) 7^ 0. 
Assume that the Killing form of the Lie algebra is such that {b\a) = implies that there 
exists c G M" such that b = a x c. Then a x u = and {u\a) =0 together imply u = 0. 



Proof 

Let d be arbitrary in R". Write 



d = + b, (23) 

{a\a) 

with {b\a) = 0. By assumption one can write b = a x c so that 

d = -l^a + a X c. (24) 
{a\a) 

Now calculate, using a x u = and {u\a) = 0, 

{d\u) = {a X c\u) = {c\u X a) = 0. (25) 
Since d is arbitrary and the Killing form is non-degenerate there follows that u = 0. 

□ 

Proposition 1 Let be given a special solution Ps{t) of the form 

p,(t) = C + ^a,(t)5„ (26) 

j 

where C commutes with all Sj. Assume that (a(0)|a(0)) 7^ and that {b\a{0)) = implies 
that there exists c G M" such that b = a{0) x c. Assume also that the Killing form is 
non- degenerate. Then one has 

H{t) = K{t) + a{t){p,{t) - C}. (27) 

with 

, , {a\h — k) 



Proof 

Note that 

a = hxa = kxa. (29) 

Let r = h — k. Then one has = r x a. Let u = r — a{t)a with a{t) given by (128|) . Then a 
short calculation yields {u\a) = and a x u = 0. By the previous lemma this implies that 
u = 0. Hence, (EZD follows. 
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□ 

Since Ps{t) and H{t) are known, the Proposition ensures that the difference between H(t) 
and K(t) is proportional to the special solution ps(t) minus a constant operator. Then, (fT6l) 
implies 

^/*(Ps(0)) = V*{t){K{t)-H{t)}V{t) 

= -a{t)V*{t){p,{t)-C}V{t) 

= -a{t){pM-C}. (30) 
Therefore, one finds that ft{u) is the linear function given by 

ft{u) = T{t)u with T{t) = - [ dsa{s). (31) 

Jo 

Hence, (fTSll becomes 

U{t) = e^"W^=(*V(t) = V{t)e'^^'^P^^°l (32) 

This completes the second step of the method. All it needs from the first step is the knowledge 
of V{t) and the value of the Killing form {k\a). 
Introduce coefficients Vmjif) defined by 

V{t)S,V*{t) = Y,VraA-^)Sm. (33) 
m 

They satisfy 

«m(t) = ^v^j{t)aj{Q). (34) 
j 

The aj{t) are known functions. We have to find Vmj(t), solving (IMIl . such that unitary 
operators V{t) exist which solve fl33l) . 
Note that 

^{a\a) = 2{h x a\a) = 2{h\a x a) = 0. (35) 

This means that {a\a) is constant in time. Hence, it is obvious to look for solutions Vkj{t) of 
f lM|) which are matrix elements of a matrix W{t) that leaves the Killing form invariant, in 
the sense that 

{Wb\Wc) = {b\c) for all 6, c G M". (36) 

For the construction of a matrix W{t) that leaves the Killing form invariant and that satisfies 
( !34|) one has to rely on the large amount of knowledge about automorphisms of Lie algebras 
(Note that automorphisms of the Lie algebra leave the Killing form invariant). 

The next step is the construction of the operator V{t). One can make use of the fact 
that the construction can be done in any faithful operator representation of the Lie algebra. 
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Indeed, let Sj be the generators in such a representation (the hat is used to distinguish the 
generators Sj from their representations 5*^) and assume that 

Psit) = Vit)mV*it) (37) 

with 

n 

p,{t) = Y,a,{t)S„ (38) 

and 

V(t) = e*'''i(*)'^ie*''^^*'''^^ ■ ■ . e*''"*^*-*'^". (39) 

Choose for instance a V{t) of the Wei- Norman form (fT2|) Then the operators V{t), 
defined by 

V(t) = e*^i(*)'^ie*''^*-*'''^^ ■ ■ ■ e*''"*^*^'^", (40) 

satisfy 

Psit) = Vit)p{0)V*it). (41) 

The argument is that, when calculating V{t)p{0)V*{t), using the Baker-Campbell-Haus- 
dorff formula, one needs to know only the commutation relations between the generators 
5*^. Hence, the expressions for V{t) and V{t) have the same form. Note that, if then the 
generators 5*-,- are self-adjoint and the coefficients rj{t) are real, then V{t) is unitary, as 
requested. 

Factorisations of V{t), other than that of Wei and Norman, will be used below as well. 
However, the above argument can be used in all cases. 

4 SU(2) 

Let be given operators 5*1, 5*2, 5*3 that satisfy the commutation relations 

[>S'i,5'2] = iSs and cychc permutations. (42) 

The non-vanishing structure constants are = 0^231 = ^^^312 = 1 and UJ132 = UJ213 = 1^321 = 
— 1. The Lie bracket is 

(a X 6)3 = 0162 — ^2^1 and cyclic permutations. (43) 

The Killing form is 

< a\b >= -2 ^ ajbj. (44) 
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A well-known identity is 



2x X (x X = {x\y)x — {x\x)y. (45) 

This implies that the conditions of the Proposition [1] are satisfied, provided that a(0) 7^ 0. 
Indeed, assume (6|a(0)) = 0. Then (H5l) implies that h = a(0) x c with 

c = -2 //°^,",V (46) 
(a(0)|a(0)) ^ ^ 

Let be given a special solution Ps(^) of the von Neumann equation. Then we have to 
construct automorphisms W(t) of the Lie algebra such that (fT3|) holds. In addition, we have 
to calculate {k\a). The rotations of are the natural candidates for the W{t). Several 
parametrisations are possible and are treated below. 

4.1 Magnus expansion 

Any rotation of M.^ can be denoted R{n, (p), meaning a rotation around an axis determined 
by the vector n, satisfying |?t,| = 1, and by an angle (j). One has 

R{n, (p)u = cos(0)m — sin(0)a x u + (1 — cos(0))(m ■ n)n. (47) 

A well-known projective representation of the group 5*0(3) in SU{2) is 

V{n, 0) = exp ^-i^ ^ j , (48) 

where the cTq are the Pauli spin matrices. It satisfies 

V{n, (j)){a ■ a)V{n, —(f)) = cos(0)a ■ a — sin(0)(a x n) ■ cr 

+ (1 — cos(0))(a ■ n)n ■ a 
= {R{n, 0)a) ■ a. (49) 

Hence, the map a — > |a • cr is an operator representation of the Lie algebra SU (2) (note that 
aj = 2Sj). 

Let n{t) and (f){t) be smooth functions such that a{t) = R{n,(f))a{0). Then (HHll is an 
explicitly constructed representation, satisfying our requirements. One concludes that 

V{t)=exp (^-i^J^n^S^ . (50) 

This ends the first step of the method. The disadvantage of the Magnus expansion is that 
the expression for the operator K is rather complicated. The expansion is therefore less 
useful and further application of the method is omitted here. 
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4.2 Wei-Norman expansion 

A rotation in can also be described as a sequence of rotations about each of the principal 
axes. Let 



/ 1 u u 

Ri{Qi) = cos(gi) sin(5i) 

y — sin(gi) cos(gi) 

(cos(g2) -sin(g2) 
1 

sin(g2) cos(g2) / 

(cos(g3) sin(g3) \ 
-sin(g3) cos(g3) . (51) 
ij 

Then one has 

(52) 

One concludes that 

T/(i) = e^*l^le^«2^^e^*3^^ (53) 

with time-dependent angles qj{t) determined by a{t) = i?i(5i)i?2(?2)-R3(?3)a(0). The Hamil- 
tonian K is found to be 

K = -qiSi - q2{cos{qi)S2 - sin(gi)S'3} 

+g3{sin(g2)'S'i + sin(gi) cos(g2)'S'2 - cos(gi) cos(g2)'S'3} 



(54) 



4.3 Euler angles 



Another well-known parametrisation of rotations is that of Euler. It describes the rotation 
of a reference frame. When describing a rotation in the first rotation can be used to bring 
the initial vector in the plane orthogonal to direction 1. Let 

a{t) = Rs{ilj)R,{e)Rs{(l))a{Q). (55) 

Then one has for arbitrary e 

eiV'^3g§9^igt</'-3(„.^)g-^</.<.3g-feaig-iv-a3 _ {Rs{^)R,{e)Rs{(l))u) ■ a. 

(56) 

One concludes that 

V{t) = e'^^'e'^^'e''!'^', (57) 
with time-dependent Euler angles determined by a{t) = i?3(-i/')i?i(6')i?3(0)a(O). 
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The Hamiltonian K is found to be 

K = - ^(cos(^)5i - sin(V^)^2) 

-0{cos(^)53 + sin(^) cos(V^)^2 + sin(^) sin(^)^i}. (58) 

Note that one can replace (157|) by the alternative expressions 

V[t) = e*('/'+<^)'S'3g«e(cos(0)Si+sm(0)S'2) ^gg^ 

or 

These expressions might be more convenient in certain applications. 
4.4 Calculation of the time-dependent angles 

The calculation of the time- dependent angles ip, 6, may involve complicated expressions. It 
is therefore noteworthy that one can redefine the Euler angles in such a way that fl571) may 
be replaced by 

V{t) = e''l''^t)S3Qimsi^-i0{o)Si^-ims3_ (gX) 
This expression involves only two angles 0(t) and 0{t), determined by 

sin(0) = — , cos(0) = — , 

sin(^^) = cos(0) = |, (62) 



where z{t) = a/ ai(t)^ + a2(t)^ and A = y z{ty + a'^{ty. To see this, note that 
i?i(— 6')i?3(— 0) rotates the arbitrary vector a into the fixed vector A(0, IjO)""". 

Using = — (ii(i2 Q _ — 3 ^j^g Hamiltonian K{t) can be calculated as follows 

z^ z 

dV 

Kit) = ^^V*{t) 

= -6{cos{(j))Si - sin(0)5'2} - 05*3 

= -^S^ - + -2 (63) 

Using a = h X a one sees that this expression is of the form k = h — a{t)a with the function 
a{t) given by 

aihi + a2h2 

a = . (64) 

The final result is then (see (1321) ) 

Uit) = Vit) exp (^trit) J2 ^AO)S^ , (65) 
with V{t) given by (16T1) and with 

r(t) = - f dsa{s) (66) 
Jo 
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4.5 General solution of the equation x = h x x 



For the applications, discussed in Part II of this paper, it is of interest that one can now 
write down the general solution of the equations x = h x x in terms of the special solution 
a and the effective time T(t). Indeed, one has 




x{t) = A \ a2 cos(Ar + C) — ai H sin(Ar + C) 0302 



(67) 



with A, B, C arbitrary constants. The verification of this results is a matter of a straightfor- 
ward calculation. 



5 SU(1,1) 

Consider operators Si, S2, S3 satisfying the commutation relations 

[Si,S2] = iSs 
['S'2,5'3] = -iSi 

[S,,Si] = -iS2. (68) 

Introduce S[ = iSi, S2 = iS2, and ^3 = —5*3. Then S[, S2, S'^ satisfy the commutation 
relations of SU(2). Hence, the results of the previous Sections can be used to solve step 1 of 
our method. However, in step 2 the assumption is made that the generators Si,S2,S3 are 
self-adjoint. Therefore, the results of Section H] have to be rederived here for SU(1,1). 

The non-vanishing structure constants are u;i23 = 0J321 = ^132 = 1 and 00213 = ^231 = 
"-1^312 = — 1- The Lie bracket for vectors in is 



The Killing form is 



a X 6)1 


= a3&2 


- a2h 


a X 6)2 


= aibs 


- ash 


a X 6)3 


= aib2 


- a2&i 



{a\b) = -2aibi - 2a2&2 + 20363. 



(69) 
(70) 



The identity fHS]) still holds. Hence, the Proposition [T] predicts the form of the relation 
between the unitary operators V{t) and Uit). 

5.1 Representation in 

Note that the cyclic permutation symmetry of SU(2) is lost. Hence, from the Wei-Norman 
expansion 0531) and the expansion based on Euler angles (15 7p one can derive two times three 
different expansions. Let us consider just one of these, corresponding with the choice of signs 
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(-,-,+) in the Killing form, as found in (ITOl) . Then the rotations around the third axis -R3(0) 
and -R3 ('?/') remain automorphisms. But Ri{9) must be replaced by 

/I 0\ 

Pi(x) = cosh(x) sinh(x) . (71) 
\ sinh(x) cosh(x) / 

Assume now that the coefficients aj{t) are written as a{t) = i?3(?/')Pi(x)/23(0)a(O). Then 
the unitary operators V{t) are given by 

Indeed, for arbitrary u in is 

g#53gie5igi053(^ ■ a)e-''^^'e-''^'e-'^^' = {Rs{-tlj)Pi{-e)Rs{-(f))u) ■ a. 

(73) 

The corresponding Hamiltonian is 

^ = '^^*it) = ^53 + x(cos(^)5i-sin(^)52) 
at 

+0{cosh(x)5'3 + sinh(x) cos{ip)S2 + sinh(x) sin (■?/;) Si}. (74) 

As before in the SU(2) case one can rewrite the automorphism in such a way that only 
two angles are involved 

ait) = P3(0(t))Pi(x(t) - xmM-m)- (75) 
The corresponding unitary operators are 

The angles (f>(t) and x(^) must satisfy 

sin(0) = — cos(0) = — 

sinh(x) = — cosh(x) = -, (77) 
/i /i 



where z{t) = y ai{ty + a2{ty and /i = y z{tY — a-i{tY. The fixed vector is 

Pi(-x(O))i?3(-0(O))a(O) = /i(0, 1, 0)^. (78) 
Using X = 0,3/ z and = (0201 — 010.2)/-^^ the Hamiltonian K can be calculated as follows 

K = t^V*{t) = x{cos(0(t))5i-sin(0(t))52} + 053 

- ^^'^i-^^'^2 + ^3- (79) 
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Using a = h X a one sees that this expression is of the form k = h — aa with the function 
a{t) given by 

aihi + a2h2 

a = . (80) 

The final result is then (see fl32l) ) 

U{t)=V{t)exp |^zr(t)^a,(0)^,j , (81) 
with V{t) given by (JTHD and with 

r(t) = - [ dsa{s) (82) 
Jo 

If the initial conditions satisfy a"^ = + a\ then = and the angle xi^) cannot be 
determined by (!77|) . In that case _R3(— 0(O))a(O) equals 2;(0, 1, 1)""" (we assume that 03(0) > 0). 
The angle x(0) can be taken equal to zero. The corresponding expression for the time- 
dependent angle xif) is then 

e^^W = (83) 

Using X = 0.3/ 0'3 and = (0201 — aia2)/al the Hamiltonian K becomes 

K = ^[a2bi - aib2) A 5 63- (84) 

One then verifies that k = h — aa still holds with the function a{t) given by (ISO]) , in spite of 
the fact that in this case the conditions of the Proposition [1] are not satisfied. 

Note that slightly different choices have to be made in the case that the initial conditions 
satisfy > of + instead of < a{ + a^. Because 03 (t) does not change sign, and is 
assumed to be positive, one can take 

sinh(x) = - cosh(x) = — , (85) 
/i /i 



with /i = Va3(t)2 - z{tf. Then Pi(-x(O))i?3(-0(O))a equals /i(0, 0, 1)^. 
5.2 General solution of the equation x = h x x 

In the same way as for the SU(2) symmetry one can now write down the general solution of 
the equation x = h x x. One finds, assuming 7^ + Og, 

uB f ^"^ \ B I "'^"'^ 

x{t) = A \ a2 I + — cosh(/ir + C) ( -ai j H sinh(/ir + C) ( 0302 | , (86) 

with A, B, C arbitrary constants. In the case that a\ = a\ + a| then one finds 

.(t) = (A + + Cr') \a.\^%(B^ 2Cr) { Z ) + | -a. | ■ (87) 
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6 Discussion 



We present a general method to construct unitary operators U (t) , solving the quantum time 
evolution in the case of a time-dependent Hamiltonian H{t). The method starts from a 
special solution of the von Neumann equation. It generalises that of [3j, where the special 
solution is called a time-dependent invariant. Here we are interested in deriving the general 
solution from the knowledge of a special solution. The method is worked out in detail for 
the case that H{t) can be written as a linear combination of generators 5*1, ■■ ■ , of a Lie 
algebra, with time-dependent coefficients. Most studied in the literature is the case of SU(2). 
The explicit expressions for the operators U{t), known as Magnus expansion, respectively 
Wei-Norman expansion, are shown to be special instances of a more general theory and 
correspond to specific representations of the rotation group in M^. A simpler expression for 
the operators U {t) is obtained when the rotations are described in terms of Euler angles - 
see ( 16T1) and (!65l) . Also the case of SU(1,1) has been discussed often. The representation of 
the rotation group involving Euler angles can be adapted to this case. The result is given 
by (Eg EH). 

The generator corresponding with the special solution is denoted K{t). The result of 
Proposition 1 is very convenient because it proves that the difference between the Hamil- 
tonians H{t) and K{t) is a linear function of the special solution Ps{t). The conditions of 
this Proposition may not always be fulfilled. But one can fall back on the general method 
described in Section [2l It suffices that the spectrum of the special solution ps(t) is discrete 
and non-degenerate. Then the unitary operators U{t) and V{t) are linked via f[T^ . The 
determination of the function ft{x) is in this case only slightly more difficult than when the 
Proposition can be applied. 

Initially [H [2], the Wei-Norman method translated the problem of obtaining explicit 
expressions for the operators U{t) into sets of differential equations with time-dependent 
coefficients. In the Lie-algebraic context these equations can be written into the form x = 
hx X, where we use the notation hx x for the Lie bracket of h and x. The main result of the 
present approach is that we obtain the general solution x of these equations, starting from 
a special solution - see fl67j) and fl86l [87|) . 

The only Lie algebras discussed in detail in the present work are SU(2) and SU(1,1). The 
extension to other finite-dimensional semi-simple Lie algebras seems to be straightforward. 

The Schrodinger and von Neumann equations with time-dependent Hamiltonian are re- 
lated to the non-linear Schrodinger and von Neumann equations. In particular, a solution 
of the non-linear equation 

^p(t) = ^[p(^)^i^] (88) 

can be used as a special solution of the von Neumann equation (j4j) with time-dependent 
Hamiltonian H{t) = p{t)H + Hp{t). See for instance [H]. The theory of finding special 
solutions for non-linear von Neumann equations can be found in jlOj . 
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